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Abstract 

We examine the feasibility of using medium resolution spectra for determining the parameters 
of atmospheres of hot stars by means of numerical simulations. We chose the star Cyg OB2 N 2 7 as 
a test object and obtained its spectrum (A/AA = 2500) with Russian- Turkish RTT150 telescope. 
The CMFGEN code was used to construct a model of the atmosphere of Cyg OB2 JY fi 7 . For the 
first time we have detected the NIV AA7103.2 — 7129.21ines in the spectrum of this star and used 
them to determine the physical parameters of the wind. The rate of mass loss measured using the 
Ha line exceeds the loss rate measured using lines formed in the wind. This indicates that the wind 
is nonuniform, apparently due to rotation. 
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1 Introduction 

Since the 1960's, space based observations made it possible to detect lines with a P Cyg profile in the 
ultraviolet spectra of O-type stars, which indicated the presence of supersonic winds. This discovery 
allowed to estimate the mass loss rate due to a stellar wind. Further studies showed that massive stars 
(> 50 Mq) lose a substantial fraction of their mass (almost half) in the form of wind while on Main 
Sequence. Later, wind lines (or details of their profiles) were also discovered in high-resolution spectra 
in the optical wavelength range obtained at ground-based observatories. 

The next important step for astrophysics was the discovery of a new method for determining dis- 
tances using the so-called wind momentum- luminosity relation (WLR) (Kudritzki et al. 1999). However 
spectral monitoring of selected O-supergiants has shown that the winds from hot stars are inhomoge- 
neous and vary with time (see Owocki (1994) for details). Therefore, to determine their luminosities it is 
important to obtain averaged, statistically significant characteristics of details of the wind line in spec- 
tra of rare O-supergiants. Significant fraction of these stars can only be studied with medium-resolution 
spectroscopy. Moreover blue parts of their spectra may be attenuated as a result of interstellar and 
circumstellar extinction. 

In this paper we test whether is it possible to reliably determine the atmospheric parameters by 
modeling medium-resolution spectra. We have chosen the supergiant 031/* Cyg OB2 N 2 7 as a test 
object. Blue part of its spectrum is significantly absorbed, A v = 5.4 (Kiminki et al. 2007). The star 
belongs to the Cyg OB2 association (Klochkova et al. 2011; Klochkova & Chentsov 2004), so we can 
obtain an independent estimate of its luminosity. 

In the next section we present observational data and data reduction process. In Section 3 we 
describe model calculations and discuss results, comparing these with previous work. The conclusions 
are considered in Section 4. 
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2 Observations and data reduction 

Observations of Cyg OB2 N 2 7 have been performed on February-March 2012 on the 1.5-m Russian- 
Turkish RTT-150 telescope at the Tiibitak National Observatory, located on Mt. Bakyrlytepe (Turkey). 
The spectrum was obtained over a wide range of wavelengths (4200-8000 AA) using the TFOSC 
(Tiibitak NationalFaint Object Spectrograph and Camera 1 ) instrument at the Cassegrain focus. The 
spectral resolution is A/ A A = 2500. In the overall spectrum the signal-to- noise ratio S/N = 100 in the 
blue part (5000 A) and 200 in the red part (7000 A). As noted above, the spectrum of this object is 
strongly absorbed at shorter wavelengths. Reliable modelling of faint lines requires a fairly high S/N 
ratio (>100) which, in our case, was realized only in the red for A > 5000 A. Below we study the 
spectrum in this region. Data were reduced using the DECH software package . 

The CIV AA5801.3, 5812, NIV AA6214, 6219 emission lines are present in the spectrum. For the 
first time the NIV AA7103.2— 7129.2 lines were detected in the spectrum of Cyg OB2 N 2 7 , they involve 
a transition from the ls 2 2s3d to the ls 2 2s3_ state. These emission lines are typical for spectra of early 
Wolf-Rayet (WR) stars and are used for a spectral classification of nitrogen-rich WR (WN). From 
published spectra of O-stars in the 7000-8000 AA range we may conclude that the NIV AA7103.2 — 
7129.2 lines are present only in the spectra of 02-05 supergiants. Modeling shows that these lines are 
formed only when effective temperature > 38000i^ (Maryeva et al. 2012). 

3 Modeling results 

We have used the CMFGEN atmosphere code (Hillier & Miller 1998) to determine the physical param- 
eters of the atmosphere of Cyg OB2 N s 7 . This code solves radiative transfer equation for objects with 
spherically symmetric extended outflows using either the Sobolev approximation or the full comoving- 
frame solution of the radiative transfer equation. CMFGEN incorporate a line blanketing, the effect 
of Auger ionization and clumping. Every model is defined by a hydrostatic stellar radius R*, lumi- 
nosity L*, mass-loss rate M, filling factor /, wind terminal velocity stellar mass M, and by the 
abundances Zi of included elementary species. 

Using the model of the star AV 83 (07 Iaf + ) calculated by Hillier et al. (2003) as the seed model 
we have adjusted its parameters to reproduce the observed spectrum of Cyg OB2 N 2 7 and gradually 
changed the parameters of the model (L*, R* and M). 

In our calculations we assumed that: 

- the volume filling factor is equal to 0.1, as in the initial model; 

- /3— law for wind velocity and Voo = 3080 km s -1 (a value taken from Herrero et al. (2001)); 

- H, He, C, N, O, S, Si, P, and Fe were included in calculations; 

- the abundances of S, Si, P, and Fe are solar abundances; 

- the abundances of H, He, C, N, and O are the same as in the initial model, ([X(N)/X(N)q] ~ 3, 
[X(C)/X(C) Q ] = 0.08, [X(0)/X(0) Q ] = 0.09). 

We used photometric data for an exact determination of the luminosity. The model flux was 
recalculated for the distance of Cyg OB2 association (1.5 kpc, Mel'nik & Dambis (2009)). Then we 
added interstellar absorption using the IDL program FM-UNRED (W. Landsman) which uses the 
absorption curves calculated by Fitzpatrick [12]. The value A v = 5.4 is taken from Kiminki et al. 
(2007). Then, the simulated spectra were convolved with V-band sensitivity filters. The resulting 
fluxes were converted to magnitudes and compared to the photometrical data (V = 10. m 5, Simbad 
data base). 



1 www. tug. tubitak.gov.tr/rttl50_tfosc.php 
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Figure 1: Left: A comparison of the observed Ha + Hell A 6560 (solid line) spectrum with the models 
reported here. The dotted line is a model with (3 = 1 and M c i = 2 x 10 _6 MQ/year and the dashed 
curve, with f3 = 1 and M c \ = 2.5 x 10~ 6 MQ/year. Another line, HellA 6527, is seen to the left of 
Ha and DIB A6613 to the right. Right: A comparison of the theoretical and observed profiles of 
Ha + HellA 6560 from Herrero et al. (2002). 



As a result, we constructed a best-fit model. Its parameters are: L* = (1.1 ± 0.1) x 10 6 Lq, R* = 
16.5i?0 , T* = 44±l/cK" and ft = 1. is the radius of the star, corresponding to the inner boundary of 
the atmosphere lying approximately at r ~ 20, and T* is the effective temperature at radius i?*, which 
is related to the luminosity by L* = AitR^aT^. The mass loss rate is M c i = (3 ± 0.5) x 10 _6 M Q /year. 
The unclumped mass loss rate (M unc i) is related to the clumped rate (M c /) by M unc i = M c ij\f] . 

Figure 1 shows comparison between the calculated and observed Ha + Hell A6560 profiles. The 
rotational velocity of the star is V smi = 105 km s" 1 (Herrero et al. 2002). In order to account for the 
star's rotation and the spectral resolution of the instrument (AA = 2A), we convolved the calculated 
spectrum with a Gaussian of FWHM=2.65A. 

In order to describe the profiles of the CIVAA5801.3, 5812 and NIV AA7103.2 - 7129.2 wind lines, 
we had to construct a model with slower acceleration of wind, corresponding to larger values of (3 
((3 = 2) and a lower mass loss rate M c \ = (3 ± 2) x 10 -7 Mo/year (Figure 2). 



Table 1: Derived properties of Cyg OB2 Wl . 
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* - The data were taken from (Herrero et al. 2002) 
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Figure 2: The profiles of NIVAA7103.2 - 7129.2 (left) and CIVAA5801.3, 5812 (right). Our model 
(P = 2, M = 6 x 10~ 7 M o /year) is shown as the dot-dashed line. DIB AA5797.03, 5809.10 appear in 
the observed spectrum. 

We now compare our results with earlier studies of Cyg OB2 N 2 7 . Herrero et al. (2002) studied its 
spectrum over a wide range of wavelengths (4000-6700 AA) using the FASTWIND code (Santolay-Rey, 
Puis & Herrero 1997; Puis 2005) which includes the blanketing effect. Their results for the Ha line are 
shown on the right panel of Figure 1. Table 1 lists the parameters of our models and the parameters 
obtained by Herrero et al. (2002). -R2/3 i s t ne radius at which the optical depth r becomes equal 
to 2/3 and T e ff is the effective temperature of the object at R2/3 (assuming radiative equilibrium). 
From the table it can be seen that our estimates for the Ha line are in good agreement with the 
earlier measurements Herrero et al. (2002). Note that the nitrogen and carbon emission lines were not 
modelled in Herrero et al. (2002), so that the differences we have found in the parameters derived from 
the Ha line and the wind lines do not contradict earlier results. 

The rotation of early-type stars with radiatively driven winds leads to interesting effects, the most 
prominent is the tendency to concentrate the outflowing material toward regions near the equatorial 
plane. This results in a deviation from a spherically symmetric shape and possibly the formation of 
an outflowing disk (Lamers & Cassinelli 1999). Disks of this kind have been discovered in B[e] and Be 
stars (Zickgraf et al. 1985). Asymmetric winds were found in objects that are evolutionary related to 
O-stars: Luminous Blue Variables (LBV) (Groh et al. 2006, 2010) and Wolf-Rayet stars (Harries et al. 
1998). Moreover numerical calculations Groh et al. (2008) have shown that density of the wind of qWR 
star HD 45166 varies with latitude. We assume that the difference between the model describing Ha 
and the one describing CIV and NIV lines is related to a latitudinal inhomogeneity of the supergiant 
wind due to its rotation (V s'mi = 105 km s _1 ). 

4 Conclusions 

We studied one of the hottest stars in our Galaxy Cyg OB2 N 2 7 by spectrum obtained on the RTT150 
telescope. Using non-LTE code CMFGEN we estimate its parameters (bolometric luminosity, stel- 
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lar radius, mass loss rate, wind velocity, elementary abundances). The atmosphere of this object 
is rich with nitrogen. We have shown that the wind of Cyg 0B2 N 2 7 is inhomogeneous. Therefore 
Cyg 0B2 N 2 7 is yet another star with density of wind depending on latitude. 

The good agreement between the parameters of Cyg OB2 N 2 7 found here and those determined from 
spectra over a wide range of wavelengths (Herrero et al. 2002) indicates that medium-resolution spectra 
in the red can be used to obtained fairly accurate estimates for the parameters of the atmospheres of 
hot stars when used in combination with reliable codes such as CMFGEN. Spectra of hot stars at red 
wavelengths contain lines formed in a stellar wind. Thus, when a red-sensitive detector is available, 
medium-resolution spectra can be used to study the strong photospheric lines, as well as features of 
the wind by monitoring its spectral variability. 
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